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Abstract 

We show that high-energy neutrinos can be efficiently produced in X-ray binaries with 
relativistic jets and high-mass primary stars. We consider a system where the star 
presents a dense equatorial wind and the jet has a small content of relativistic protons. 
In this scenario, neutrinos and correlated gamma-rays result from pp interactions and 
the subsequent pion decays. As a particular example we consider the microquasar LS I 
+61 303. Above 1 TeV, we obtain a mean-orbital z^-luminosity of ~ 5 10 34 erg/s which 
can be related to an event rate of 4-5 muon-type neutrinos per kilometer-squared per 
year after considering the signal attenuation due to maximal neutrino oscillations. The 
maximal neutrino energies here considered will range between 20 and 85 TeV along the 
orbit. The local infrared photon-field is responsible for opacity effects on the associated 
gamma radiation at high energies, but below 50 GeV the source could be detected by 
MAGIC telescope. GLAST observations at E~ > 100 MeV should also reveal a strong 
source. 



PACS: 95.85.Ry; 98.70.Sa; 13.85.Tp 13.85-t 

To Appear in Phys. Rev. D (April 2006 issue) 1 



1 Introduction 

The study of high-energy neutrinos from galactic sources is expected to provide important 
clues for the understanding of the origin of cosmic rays in our Galaxy. Astrophysical sources 
of high-energy neutrinos should have relativistic hadrons and suitable targets for them, 
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such as radiation and/or matter fields. Rotating magnetized neutron stars are well-known 
particle accelerators. When they make part of a binary system, relativistic particles may 
find convenient targets in the companion star and high-energy interactions can take place. 
It is then natural to consider accreting neutron stars or, more generally, X-ray binaries, as 
potential neutrino sources (see Bednarek et al. 2005 and references therein^)- 

Neutron stars with high magnetic fields can disrupt the accreting flow, which is channeled 
through the field lines to the magnetic poles of the system. The magnetosphere of such 
systems presents electrostatic gaps where protons can be accelerated up to very high energies. 
Actually, even in the absence of relativistic jets, it is possible to attain hundreds of TeVs 
(Cheng & Ruderman 1989, 1991 [21 IS])- The accelerated protons, which move along the 
closed field lines, can impact on the accreting material producing gamma-rays (Cheng et al. 
1992a [1], Romero et al. 2001 Orellana & Romero 2005 [B]) and neutrinos (Cheng et al. 
1992b [7], Anchordoqui et al. 2003 0). 

If the magnetic field of the neutron star is not very high, then accretion/ejection phenom- 
ena can appear, and the X-ray binary can display relativistic jets, as in the well-known cases 
of Sco X-l (Fender et al. 2004 [SJ) and LS I +61 303 (Massi et al. 2001, 2004 [HI EE]). The 
X-ray binary is then called a 'microquasar'. In this case, a fraction of the jet hadrons can 
reach much higher energies, up to a hundred of PeVs or more, depending on the parameters 
of the system (see below). 

Microquasars can be powered either by a weakly magnetized neutron star or by a black 
hole (as, for instance, Cygnus X-l). Some of them are suspected to be gamma-ray sources 
(Paredes et al. 2000 [121, Kaufman-Bernado et al. 2002 03], Bosch- Ramon et al. 2005a [13]). 
Recently, Aharonian et al. (2005) J3] have detected very high-energy emission from the 
microquasar LS 5039 using the High Energy Stereoscopic System (HESS). In microquasars, 
the presence of relativistic hadrons in the jets can lead to neutrino production through 
photo-hadron (Levinson & Waxman 2001 [TE] . Distefano et al. 2002 [17] ) or proton-proton 
(Romero et al. 2003 HE], Romero & Orellana 2005 Torres et al. 2005 [20], Bednarek 
2005 [21]) interactions. In the latter case, the target protons are provided by the stellar wind 
of the companion. 

In the case studied by Romero et al. (2003) [IE], which was the base of subsequent models, 
the wind is assumed to be spherically symmetric. This, however, is not always the case. In 
particular, microquasars with Be stellar companions present slow and dense equatorial winds 
that form a circumstellar disk around the primary star and the compact object moves inside. 
In the present paper we will study how the interaction of the jet with this material can lead to 
a prominent neutrino source. We shall focus on a specific object, LS I +61 303, for which all 
basic parameters are rather well-determined in order to make quantitative predictions that 
can be tested with the new generation of neutrino telescopes (IceCube, in this case |22j). 
Our results, however, will have general interest for any microquasar with non-spherically 
symmetric winds. We shall start with a brief description of LS I +61 303. 

2 The microquasar LS I +61 303 

LS I +61 303 is a Be/X-ray binary system that presents unusually strong and variable radio 
emission (Gregory & Taylor 1978 |23j). The X-ray emission is weaker than in other objects 
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of the same class (e.g. Greiner & Rau 2001 [21]) and shows a modulation with the radio 
period (Paredes et al. 1997 [22] )• The most recent determination of the orbital parameters 
(Casares et al. 2005 j2H]) indicates that the eccentricity of the system is 0.72 ± 0.15 and 
that the orbital inclination is ~ 30° ± 20°. The best determination of the orbital period 
(P = 26.4960 ± 0.0028) comes from radio data (Gregory 2002 [2B])- The primary star is a 
B0 V with a dense equatorial wind. Its distance is ~ 2 kpc. The X-ray/radio outbursts are 
triggered 2.5-4 days after the periastron passage of the compact object, usually thought to 
be a neutron star. These outbursts can last until well beyond the apastron passage. 

Recently, Massi et al. (2001) ^Oj have detected the existence of relativistic radio jets 
in LS I +61 303. These jets seem to extend up to about 400 AU from the compact object 
(Massi et al. 2004 |2S!). 

LS I +61 303 has long been associated with a gamma-ray source. First with the COS- 
B source CG135+01, and later on with 3EG J0241+6103 (Gregory & Taylor 1978 [23], 
Kniffen et al 1997 j30j). The gamma-ray emission is clearly variable (Tavani et al. 1998 |3*T] ) 
and has been recently shown that the peak of the gamma-ray lightcurve is consistent with 
the periastron passage (Massi 2004 JI]), contrary to what happens with the radio/X-ray 
emission, which peaks after the passage. A complete an updated summary of the source is 
given by Massi 2005 (221 

The matter content of microquasar jets is unknown, although in the case of SS 433 iron 
X-ray line observations have proved the presence of ions in the jets (Kotani et al. 1994, 
1996 [HBIE1!; Migliari et al. 2002 [SS])- In the present paper we will assume that relativistic 
protons are part of the content of the observed jets in LS I +61 303. In the next section 
we will describe the basic features of the model (see Romero et al. 2005 [SH] for a detailed 
discussion of the gamma-ray emission), and then we will present the calculations and results. 



3 Model 

The system under consideration consists of two objects. The primary is a B-type star that 
generates a radially outflowing wind. The other is a compact object (CO hereafter, probably 
a neutron star) moving around in a Keplerian orbit. Their relative position is given by 
r(i/j) = a(l — e 2 )/(l — e cos(ip)), where ip is the orbital phase, a is the semi-major axis of the 
ellipse and e its eccentricity 2 . The primary presents a nearly equatorial circumstellar disk 
with a half-opening angle <fi = 15°. Its density is given by p w (r) = p ( r /-R*)~ n an d the wind 
velocity reads v w = Vq(t / ' R*) n ~ 2 by continuity. We will adopt n = 3.2 following Marti and 
Paredes 1995 [37] (see also Gregory & Neish, 2002 |35j) . The wind accretion rate onto the 
compact object is given by 

M c= MGMf Pw (r) 

where M c is the CO mass and t> re i is its velocity relative to the circumstellar wind. The 
kinetic jet power Qj is coupled to M c by Q] = q^M c c 2 (Falcke & Biermann 1995 39 J. In 
the jet-disk symbiosis model q$ ~ 0.1 for microquasars in the low-hard state (Fender 2001 

2 Notice that here t/j is not the radio phase, which amounts 0.23 2tt at the periastron (Casares et al. 2005) 

m 
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Figure 1: A sketch of the binary system with its elements in detail. 



). Only a small fraction of the jet particles are highly relativistic hadrons (~ 10 ), and 
these are confined by the pressure of the cold particles expanding laterally at the local sound 
speed (see Bosch- Ramon et al. 2005b jH] for a detailed discussion). In the present approach, 
most of the jet power will consist of cold protons ejected with a macroscopic Lorentz factor 
r ~ 1.25 (Massi et al. 2001 |Ej). 

The jet will be a cone with a radius Rj(z) = z(Ro/z ), where z = 10 7 cm is the injection 
point and i? = z /10 is the initial radius of the jet (see Romero et al. 2003 18] and Bosch- 
Ramon et al. 2005a for additional details). The jet axis, z, will be taken normal to the 
orbital plane. In Figure ffj we show a sketch of the general situation. 

We will use a power-law relativistic proton spectrum N' p (E' p ) = K p E'~ a (the prime 
refers to the jet frame). The corresponding relativistic proton flux is J' p (E' p ) = (c/4iT)N p (E p ) 
evolving with z as J' p {E' p ) = ^ a (z / z) 2 E f p ~ a , (conservation of the number of particles is 
assumed, see Ghisellini et al. 1985 021)- Using relativistic invariants, it can be shown that 
the proton flux in the lab (observer) frame becomes (e.g. Purmohammad & Samimi 2001 



J p (E p ,9) = — {-) 



sin 2 fl + r 2 cos9- 



1/2 ' 



(2) 



The angle subtended by the proton velocity direction and the jet axis will be roughly the 
same as that of the emerging photon (9 ~ o bs)- Pb is the bulk velocity in units of c, and 
a = 2.2 characterizes the power law spectrum (see the list of parameters in Table 1). We 
adopt a value of 2.2 for the proton power-law index instead of the canonical one of 2.0 given 
by first order diffusive shock acceleration. This is in order to match the GeV gamma-ray 
spectrum observed by EGRET jH]. Deviation from the canonical value are likely due to 
nonlinear effects (see, e.g., Malkov Sz Drury 2001 jH]). The normalization constant K and 
the number density no' of particles flowing in the jet at Rq can be determined as in Romero 
et al. (2003) [IE]. 
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As long as the particle gyro-radius is smaller than the radius of the jet, the matter from 
the wind can penetrate the jet diffusing into it 3 . However, some effects, like shock formation 
on the boundary layers, could prevent some particles from entering into the jet. Actually, 
the problem of matter exchange through the boundary layers of a relativistic jet is a difficult 
one. Contrary to cases studied in the literature, where the external medium has no velocity 
towards the jet (e.g. Ostrowski 1998 [IE]), here the wind impacts on the jet side with a 
velocity that can reach significant values (hundreds of km s -1 ). A detailed study of the 
penetration of the wind into the relativistic flow is beyond the scope of the present paper, 
so we will treat the problem in a phenomenological way. This can be done by means of a 
"penetration factor" f p that takes into account particle rejection from the boundary. We 
will adopt f p ~ 0.1, in order to reproduce the observed gamma-ray flux at GeV energies, 
where opacity effects due to pair creation are unimportant. 

We will make all calculations in the lab frame, where the cross sections for proton interac- 
tions have suitable parameterizations. Neutral pion decay after high energy proton collisions 
is a natural channel for high energy gamma-ray production. The differential gamma-ray 
emissivity from 7r°-decays can be expressed as (e.g. Aharonian & Atoyan 1996 47J): 

n j/0 3 ) 

g 7 (£ 7 , 6) = A7r VA a pp (E p )^^ J p (E y , 9), (3) 

a 

(in ph s -1 sr -1 erg -1 ), where Z^l n0 is the so-called spectrum- weighted moment of the in- 
clusive cross-section and it is related to the fraction of kinetic proton energy transferred to 
the pions (see Gaisser 1990 The parameter t/a takes into account the contribution 

from different nuclei in the wind. For standard composition of cosmic rays and interstellar 
medium t]a ~ 1.4. The proton flux distribution J p (E y ) (see eq.©) is evaluated at E = E 7 , 
and <Tpp(Ep) ~ 30 x [0.95 + 0.06 log (E p /GeV)] (mb) 47J is the cross section for inelastic pp 
interactions at energy E p m 6£ n oE 1 /K, for E 1 > 1 GeV. Here K is the inelasticity coefficient, 
and £^0 represents the pion multiplicity. For 6 we will adopt a viewing angle of 9 Q ^ S = 30° in 
accordance with the average value given by Casares et al. (2005) [2*o] . 
The spectral gamma ray intensity (in ph s _1 erg -1 ) is 

I 7 (E 7 ,0)= [ d 3 f n(f)q 7 (E^0), (4) 
Jv 

where V is the interaction volume between the jet and the circumstellar disk. The parti- 
cle density of the wind that penetrates the jet is n(r) « f p p w (r)/m p , and the generated 
luminosity in a given energy band results 

L 1 (E a ^\6)= I ^ dE 1 E^ Iy(E y , 9) (5) 

Since we are interested in predicting the total gamma ray luminosity measurable on 
Earth, we must evaluate this integral above a threshold of detection which shall depend on 

3 This imposes a constraint onto the value of the magnetic field in the jet: i?j e t > E^' wmd /(ei?), where 
E%' Wlnd = m p «rci/2 is the kinetic energy of the cold protons in the slow stellar wind (for ££' wmd maximum, 
at periastron, results Bj et > 2.8 10~ 6 G). This condition, as we will see, is assured. 
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the telescope characteristics. On the other hand, we should keep in mind that E 1 cannot 
exceed the maximum gamma ray energy available from the hadronic processes described so 
far. We will address this issue in more detail in the following section. 



4 Hadronic production of neutrinos and gamma rays 

In order to discuss the hadronic origin of high energy gamma-rays and neutrinos we have 
to analyze in more detail pion luminosity and decay. Regarding gamma-rays production, 
since the branching ratio of charged pions into photons (plus leptons) is about 6 orders of 
magnitude smaller than that of 7r°'s into photons (alone) |2Z| we shall first concentrate on 
these neutral parents. Actually, high-energy neutral pions can be produced by relativistic 
proton collisions in either pp or p7 interactions depending on the relative target density of 
photons and protons in the source region (where the protons are accelerated) and on their 
relative branching ratios. On the other hand, charged pions will be responsible for high 
energy neutrinos. In any case, note that both photon and neutrino production at the source 
are closely related. The v^+v^ differential neutrino flux (dN v /dE v ) produced by the decay of 
charged pions can be actually derived from the differential 7-ray flux (dN^/dE^). Following 
Alvarez- Muniz and Halzen (2002) 48 3 , we will find the neutrino intensity and spectral flux 
by means of an identity related to the conservation of energy (see also Stecker 1979 |I9| 150]) 

max „graax 

T E^dE, = D K-^dE v . (6) 

7 1 J -'v 

Here -E™ m (E™ ax ) is the minimum (maximum) energy of the photons that have a hadronic 
origin and E™ m and E™ x are the corresponding minimum and maximum energies of the 
neutrinos. The relationship between both integrals, as given by D, depends on the energy 
distribution among the particles resulting from the inelastic collision. As a consequence, 
it will depend on whether the 7r°'s are of pp or origin. Its value can be obtained from 
routine particle physics calculations and some kinematic assumptions. Let us admit that 
in pp interactions 1/3 of the proton energy goes into each pion flavor on average. We can 
further assume that in the (roughly 99.9%) pion decay chains, 

tt + -> /i + z/ M -> e + u e u fl u fl vr~ -> prv^ -> e~z/ e z/ M z/ M tt° -> 77, 

two muon-neutrinos (and two muon-antineutrinos) are produced in the charged channel with 
energy E n /4, for every photon with energy E v /2 in the neutral channel. Therefore the energy 
in neutrinos matches the energy in photons and D = 1. 

The relevant parameters to relate the neutrino flux to the 7-ray flux are the maximum and 
minimum energies of the produced photons and neutrinos, appearing as integration limits in 
Eq. (jUJ) • The maximum neutrino energy is fixed by the maximum energy of the accelerated 
protons (i?™ ax ) which can be conservatively obtained from the maximum observed 7-ray 
energy i?™ ax . Following our previous assumptions 
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for the pp case. The minimum gamma and neutrino energies are fixed by the threshold for 
pion production. For the pp case 

(2m v + mJ 2 - 2ml / x 

771mm P " J p /q\ 

E *> =r 2m~ p ' (8) 

where T is the Lorentz factor of the accelerator relative to the observer. The average mini- 
mum neutrino energy is obtained from E™ m using the same relations of Eq. (|7J) . 

In p7 interactions we can assume that neutrinos are predominantly produced via the 
A-resonance. In 1/3 of the interactions a 7r + is produced which decays into two neutrinos 
of energy E^/4, and in the other 2/3 of the interactions a 7r° is produced which decays into 
two photons of energy E n /2. Therefore D = 4. The minimum proton energy is given by 
the threshold for production of pions, and the maximum neutrino energy is about 5% of the 
maximum energy to which protons are accelerated, which is much more than that obtained 
from pp collisions. However, when the production of neutrinos in pry collisions occurs in 
the acceleration region, the efficiency of conversion into relativistic hadrons is lower than 
necessary. This is due to the fact that the threshold for electron-positron production is 
about two orders of magnitude below that for pion production along with the fact that, 
in this case, the e ± production cross section is significantly larger. As a result, most of the 
energy from the acceleration mechanism is transferred to leptons, radiating plenty of photons 
but no neutrinos. 

On the other hand, when purely hadronic collisions are considered, the situation is 
reverted and pion production becomes the natural production channel. Relativistic pro- 
tons in the jet will interact with target protons in the wind through the reaction channel 
p + p — > p + p + ^Trovr + £ w ±(tt + + ir~). Then pion decay chains lead to gamma-ray and 
neutrino emission. Isospin symmetry, which is in agreement with Fermi's original theory of 
pion production (where a thermal equilibrium of the resulting pion cloud is assumed) relates 
the three multiplicities with an equal sign, thus we simply write £„■. 

We can go a step further and consider energy dependent multiplicities. In this case, the 
relation between energy maxima is quite different from Eq.l jTj l- According to Ginzburg Sz 
Syrovatskii (1964) [HT], for inelastic pp interactions we can obtain the gamma ray energy 
from the proton energy by 

E p « 6 UE P ) Ej/K. (9) 

The inelasticity coefficient is K ~ 0.5 since on average a leading nucleon and a pion 
cloud leave the interaction fireball each carrying half of the total incident energy 4 . For the 
energy dependent pion multiplicity we will follow the prescription adopted by Mannheim 
and Schlickeiser (1994) (33] 

& ~ (E p /GeV - l.22 f\ (10) 

Note however that for a proton energy of 1 PeV, the relation between maxima given by eq. 
(J2J would become i?™ ax ~ 380£™ ax , and accordingly E™ x ~ ^£™ ax . This relation implies 
that the 1 TeV threshold detector IceCube should actually measure some neutrinos provided 
the microquasar could accelerate protons up to the PeV range. As a matter of fact, Eq. (jlOj) 
overestimates pion luminosity at energies higher than 10 GeV and the 1/4 root should be 

4 See ref. for a discussion on this coefficient. 
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relaxed for the proton energies involved in our calculation . Since we have no data about 
its energy dependence we shall adopt a softer root growing law 5 for £ n (E p ) which leads to 
a more reliable relationship between neutrino and proton energies. For the highest proton 
energies, as given by Eq. (fT4*|) . one gets ~ j^E™* at periastron (see FigJ2b). We will 
promptly see that one can in fact expect for a significant signal to noise relationship in less 
than one year of operational time of a 'km 3 ' telescope. 

The point now is whether a microquasar might accelerate protons to such high energies. 
As mentioned in Bednarek et al. (2005) pQ, the maximum energy at which a particle of 
charge Z can be accelerated, can be derived from the simple argument that the Larmor 
radius of the particle should be smaller than the size of the acceleration region (Hillas 1984 
[56 ). If energy losses inside sources are neglected 6 , this maximum energy E^ (in units of 
10 15 eV) is related to the strength of the magnetic field B (in units of Gauss) and the size 
of the accelerating region R (in units of cm) by the following relationship 

^15 ~ 0-3 10~ 12 /3 Z Bq Rem (11) 

where (3 is the velocity of the shock wave or the acceleration mechanism efficiency Hence, 
the maximum energy up to which particles can be accelerated depends on the BxR product. 
Diffusive particle acceleration through shocks may occur in many candidate sites with sizes 
ranging from kilometers to megaparsecs. In the case of LS I +61 303, we adopt quasi-parallel 
shocks, R ~ Rj (jet radius), and B given by equipartition with the cold, confining plasma 
(see Bosch- Ramon et al. 2005b jlTj). 

For a complete calculation of the magnetic field, we will assume a magneto-hydrodynamic 
mechanism for the ejection where both matter and field follow adiabatic evolution when 
moving along the jet. In this conditions it is reasonable to adopt equipartition between the 
magnetic field energy and the kinetic energy of the jet. This leads directly to the following 
expression for the magnetic field in the jet reference frame at different distances z from the 
compact object [ST] [4*T] 



B(z) = v/87re jet , (12) 

where 

e JC t = p2 7 J f E k (13) 
7iR*(z)Vjm p 



5 According to Begelman et al.(1990) based on Orth and Buffington (1976) [S3], the multiplicity 
should be kept below 15 for E p = 100 TeV. In order to analytically lessen Eq. 110(1 to such values, we 
can interpolate the 10 to 100 TeV range in a simple way with a 1/5 fractional power. This gives £ w w 14 
for E p = 100 TeV. For the next two decades, there are no confident approaches to our knowledge so we 
extrapolate £„• with a 1/6 root up to the PeV maximal proton energies obtained in Eci. (|14J) . This gives £„• « 
26 for E p = 50 PeV (periastron). Notice that a continuous expression for £ ff is needed in order to evaluate 
both luminosity and neutrino signal from Eq. J5J) and Eq. <|19f) . Notwithstanding, the corresponding integrals 
are not really sensitive to this analytical choice since the multiplicity effects are actually weak: on the one 
hand, a pp depends only logarithmically on £„.; on the other, it modifies the upper limit of the integrals, but 
then again as the integrand is quite steep the contribution of the queue in a larger domain is below a few 
percent. Note finally that here we conservatively adopted a log instead of a log 2 Ep - dependence of the 
cross section. 

6 Somcthing that, given the size constraints and photon fields, is the case in the jet of a source like LS I 
+61 303. 
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Figure 2: Left (a): Maximal magnetic field in the jet, B max (ip), as a function of the orbital 
position of the compact object (in Gauss). Right (b): Highest proton energy compatible 
with the magnetic field (in PeV) as a function of if). 



is the jet energy density for a cold proton dominated jet. The mean cold proton kinetic 
energy, Ek, is taken to be the classical kinetic energy of protons with velocity uj = The 
ejected matter amounts to a fraction of the ratio of accreted matter, namely mj et = q^M c 
with ~ 0.1. 

Now, we can proceed and compute the maximum proton energy compatible with this 
magnetic field by means of the gyro-radius identity 7 

E™ x = R(z) eB(z). (14) 

This result will be used to compute the maximum jet gamma-ray and neutrino energies 
which we will need in the next section to calculate both intensities. 

In Figure 121 we show the magnetic field (a) and maximum proton energy (b) at the base 
of the jet in terms of the orbital position for LS I +61 303. 

As we can see in the figure, the maximum proton energy produced at the base of injection 
amounts to more than 50 PeV at periastron and lies below 15 PeV near apastron. Note that 
the minimum does not occur exactly at apastron due to the magnetic field dependence. This 
is a consequence of the particular orbital variation of the accretion rate. Furthermore, we 
see that there is a second sharp maximum, readily identified with the function l/ff el , which 
grows dramatically when the compact object gets far from the primary star and moves more 
parallel to the weak stellar wind (see also Marti & Paredes 1995 jH7] for a discussion of this 
effect). 



5 Flux and event rate of neutrinos 

In order to obtain the spectral flux of neutrinos, we need an explicit expression for the 
neutrino emissivity with the details therein. Instead, we can take eq. © so as to extract an 
expression for I V (E U ). Let us write eqs. © and (JJJ) as 

L U = D L 7 , (15) 

7 Notice that a straightforward calculation of the proton losses shows that they do not diminish the 
maximum energy since the size constraints are more important in the present context. 
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and 



2E V = E 7 , (16) 

which result from the assumption that, together, one (//) neutrino and one anti-neutrino 
carry half one charged-pion energy (note that the number of (/i) neutrinos resulting from tt + 
and 7T~ decays is equal to the number of gamma rays coming from tt° since the number of 
pions produced is the same for the three flavors; same holds for the flux of (//) anti-neutrinos). 
From the last two equations we obtain 

I V {E V ) = AD I 7 (£ 7 = 2E V ). (17) 

For pp interactions 1/3 of the proton energy goes to each pion flavor and we have D = 1. 
Thus we have for the neutrino intensity 

I U (E U , 6) = 4 / dV fp Pw{r) g 7 (f, 2E U , 9) (18) 
Jv m P 

Now, we have to compute the convolution of the neutrino flux with the event probability 
For signal and noise above 1 TeV we obtain respectively 



1 l /-27T rEfW 

S = T ohs A cS -—— dip I dE v I v {E v ,tP) P V ^{E V ), (19) 

4ird z 2tt J J 1TeY 

and 



N 



1 /-27T rE™**(i>) 

T ohs A eS AQ— / dtl) / dE v F B {E V ) P U ^(E V 

^ JO JlTeV 



1/2 



(20) 



where T Q b s is the observational time period, A e g is the effective area of the detector, d is the 
distance to the system, and AQ is the solid angle of the search bin. The function 

F B (E„) < 0.2 (£7„/GeV)- 3 - 21 GeV^cmVsr" 1 (21) 

represents the + atmospheric flux (see Volkova 1980, and Lipari 1993 [SHI EH]), an d 

P V ^{E V ) = 1.3 10- 6 (E./TeV) - 8 (22) 

is the probability that a neutrino of energy E v ~ 1 — 10 3 TeV, on a trajectory through the 
detector produces a muon (see Gaisser et al. 1995 |5D*]). 

As we have already discussed, the value of the maximal neutrino energy is a function of 
orbital position of the jet and is given by 

« T^T) ^ aX (^)- (23) 

Our estimate for this expression is £"™ ax w 85 TeV at periastron, and about 20 TeV at 
apastron (see FigEJ). 

If we now set the features of a km-scale detector such as IceCube [22] in the model 
(A e s = 1 km 2 and AQ w 3 10 _4 sr), and the distance to LS I +61 303 (d — 2 kpc) we get 

S/V| lycar > 9.42 (24) 
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1 2 3 4 5 6 

Figure 3: Maximal neutrino energies i?™ ax along the orbit (in TeV). 

for the signal to noise ratio in a one-year period of observation. This is consistent with the 
upper limit recently reported by Ackermann et al. (2005) from AMANDA-II experiment 

Note that the present prediction is restricted to /z-neutrino production at the source. As 
a matter of fact, there is nowadays strong experimental evidence of the existence of neutrino 
oscillations which occur if neutrinos are massive and mixed. In 1998 (see Super-Kamiokande 
[62 ) and in 2002 (see K2K [03]) it has been observed, respectively, the disappearance of 
atmospheric and laboratory muon-neutrinos as expected from flavor oscillations. Also in 
2002, the SNO experiment [S3] provided solid evidence of solar electron-neutrino oscillations 
to other flavors. Solar and atmospheric neutrino flux suppression can be explained in the 
minimal framework of three- neutrino mixing in which the active flavor neutrinos z/ e , and 
v T are unitary linear combinations of three mass (Majorana) eigenstates of the neutrino 
lagrangian [53] . 

The expected ratios at sources of high-energy neutrino fluxes from pp collisions are 1 : 
2 : 10~ 5 for the e, /i, r neutrino flavors, in the range 1 GeV < E u < 10 12 GeV. The neutrino 
oscillation effects imply that one should measure different values depending on the distance 
to the astrophysical source. The estimate is that these ratios become in average 1:1:1 for 
L(pc)/E(GeV) > 10~ 10 , where the distance is in parsecs (see Athar et al. 2005 jUS])- Also 
recently, a detailed analysis of supernova remnants reported by Costantini and Vissani (2005) 
[6*7] claims a 50% of muon-neutrino plus muon anti-neutrino flux reduction due to flavor 
oscillations along astrophysical distances. Note finally, that the neutrino signal measured at 
Earth might be further attenuated due to matter absorption 8 (see both references above) 9 . 

6 Luminosity and opacity 

As a result of the discussion above, our predictions for both neutrino and gamma ray lumi- 
nosities between -E 1 ™ 1 " = 1,2 TeV and -E™ a * amount to 5 10 34 erg s _1 for a source such as LS 
I +61 303. 

8 Separate (in matter) flavor oscillations effects for its path across the Earth should be also considered for 
an exhaustive analysis. 

9 According to the central values of the mixing angles reported in Costantini and Vissani 2005 it is 
not clear if there could also be a slight contribution to muon-neutrino detection due to electron-neutrino 
oscillation. As for tau-neutrinos, they should be already highly suppressed at the source. 
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The total neutrino luminosity is related to the spectral neutrino intensity by means of 



dE V E U lyiEy) 



(25) 



as we mentioned in Sect. EJ Since we do not have an expression for I U (E U ), we can use eq. 
(jHl) and relate it to the gamma ray luminosity 



^max 



(1E^ Ery IryiEr, 



(26) 



which we know in more detail. Now, there is in fact an experimental constraint on this 
quantity, as recently analyzed by Fegan et al. (2005) jHHl- There, it is claimed that, above 
0.35 TeV, the total flux of gamma rays that can be produced must satisfy 



And 2 



dE 1 I 7 (£ 7 ) < 1.7 10~ n ph cm-V 1 . (27) 

'0.35TeV 

Our calculations for gamma rays give about four times this value, so we need to show that 
most of the high-energy gamma rays get absorbed. Of course, the absorption mechanism 
should not modify the flux of neutrinos. Infrared photon fields can be responsible for TeV 
photon absorption in the source. To show it, we calculate the optical depth r within the 



circumstellar disk for a photon with energy E^ 

poo poo 

r(p, EJ = / dE ph / dp n ph (E ph , p) a e - e +(E ph , E^, 



(28) 



where E^ is the energy of the ambient photons, n p h(-Eph, p) is their density at a distance p 
from the neutron star, and <7 e - e +(-E p h, £? 7 ) is the photon-photon pair creation cross section 
given by 

,2 



&e+e-{E p h, E^ 



TIT, 



2£{e - 2) + (3 - £ 4 ) In 



1 + g 
1-^ 



where tq is the classical radius of the electron and 

(m e c 2 ) 2 



E p h_Ej 



1/2 



(29) 



(30) 



In Eq. 



j E m \ n 



m e c ./E~ is the threshold energy for pair creation in the ambient photon 



field. This field can be considered as formed by two components, one from the Be star and 
the other from the cold circumstellar disk 
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n ph:1 (E ph ,p) 



7T 



B 



(£ P h) \ 



™ P h,i + ™ P h,2- Here, 
Rl 



(32) 



hcE p h J p 2 + r 2 — 2pr sin 6* ' 

10 There is, in fact, a third component corresponding to the emission from the heated matter (hot accreting 
matter impacting onto the neutron star) which can be approximated by a Bremsstrahlung spectrum 



lph,3(^ph, p) — 



iT 2 



A-ncp 2 e^ph/^cut-ott 



for £ ph > 1 keV, 



(31) 



where Lx is the total luminosity. This component does not significantly affect the propagation of TeV 
gamma-rays so we neglect this contribution. 
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Figure 4: Optical depth at periastron, as a function of log (.Ey/ GeV). 
is the black body emission from the star, and 

f nB(E ph ) \ r 2 

n ^ E ^ p) = {-h^)7' (33) 

corresponds to the emission of the circumstellar disk. In both cases we adopt 

B{Eph) = {hcf (^W _ l) ( 34 ) 

where T eff)1 = 22500 K and T cffj2 = 17500 K (Marti & Paredes 1995 [37]). 

The mentioned value for T e fj j2 is valid in the inner region of the disk and thus reliable only 
near periastron (as one gets far from the central star the temperature of the disk gets lower). 
Hence we adopt r = r per i astr . In any case, computing the opacity at periastron, where the 
luminosity is maximal, will be enough for our purposes of proving that the total gamma-ray 
luminosity predicted in our model is below the Fegan et al.'s constraint jHS]- Figure |U shows 
the Ery dependence of the optical depth at periastron for an observer at # b s = 30° with 
respect to the jet axis. The optical depth remains above unity for a wide range of photon 
energies but a sharp knee takes place at E 1 ~ 250 GeV. For TeV gamma energies it still has 
an important effect on the luminosity. In particular, the opacity-corrected total flux above 
Ery = 350 GeV drops from 7.16 10~ n to 1.42 10 -12 ph cm~ 2 s _1 , lying well below the Fegan 
et al.'s threshold. 

Other observing windows, however, could reveal the presence of a hadronic gamma-ray 
source at the position of LS I +61 303. In particular, between 1 and 50 GeV the opacity 
is sufficiently low as to allow a relatively easy detection by instruments like the ground- 
based MAGIC telescope and the LAT instrument of GLAST satellite. The source is in the 
Northern Hemisphere, hence out of the reach of HESS telescopes, which recently detected 
the microquasar LS 5039 11 . LS I +61303 is an outstanding candidate to corroborate that 
high-energy emission is a common property of microquasars. Its location, in addition, makes 
of it an ideal candidate for IceCube. A neutrino detection from this source would be a major 
achievement, which would finally solve the old question on whether relativistic protons are 
part of the matter content of the energetic outflows presented by accreting compact objects. 

n Late after the completion of this paper we became aware of a subsequent upload 69, where the potential 
TeV neutrino source LS 5039 is considered. 
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7 Conclusions 



We have analyzed the possible origin of high-energy neutrinos and gamma rays coming 
from a galactic source. Our specific subject for the analysis has been the microquasar LS 
I +61 303 because it is a well-studied object of unique characteristics and an appropriate 
candidate to make a reliable prediction of the neutrino flux from a galactic point source. We 
performed our calculations within a purely hadronic framework and showed how neutrino 
observatories like IceCube can establish whether TeV-gamma rays emitted by microquasars 
are the decay products of neutral pions. Such pions are produced in hadronic jet - wind 
particle interactions. We improved previous predictions by considering realistic values for 
the parameters of the system and energy dependent pion multiplicities particularly significant 
at high energies. Above 1 TeV, we obtained a mean-orbital //-luminosity of 5 10 34 erg/s which 
can be related to an event rate of 4-5 muon-type neutrinos per kilometer-squared per year 
if we take into account neutrino oscillations. The upper limit of integration depends on 
the orbital position and is a function of the largest magnetic gyro-radius compatible with 
the jet dimensions. As a consequence, the maximal neutrino energies here considered range 
between 20 and 85 TeV along the orbit. Opacity effects on the associated gamma radiation 
are due to the infrared photon field at the source and result in a significant attenuation 
of the original 7-ray signal. Nonetheless, LS I +61 303 might be detectable at low GeV 
energies with instruments like MAGIC and GLAST. Such a detection, would be crucial to 
test current ideas of particle acceleration in compact objects. 
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Table 1: Basic parameters assumed for the model 


Parameter 


Symbol 


Value 


Mass of the compact object 
Jet's injection point 1 


M c 


1.4 M 




50 R g 


Initial radius 


Rq 


zo/10 


Mass of the companion star 


M* 


10 M 


Radius of the companion star 


R* 


10 R Q 


Effective temperature of the star 




22500 K 


Density of the wind at the base 


Po 


10~ n gr cm" 3 


Initial wind velocity 2 


v 


5 km s _1 


Jet's Lorentz factor 


r 


1.25 


Minimum proton energy 


771/min 


1 GeV 


Penetration factor 


U 


0.1 


Orbital axis (cube) 


a 3 


P 2 G{M i ,+M c ) 
4tt 2 


Eccentricity 


e 


0.72 


Orbital period 


P 


26.496 d 


Relativistic jet power coupling 


Vj 


0.01 


Index of the jet proton distribution 


a 


2.2 



1 /.', = GM c /c 2 . 

2 The wind velocity increases from this value up to the escape 
velocity according to the velocity law given in the text. 
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